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Abstract. The dynamics of confined glassforming liquids is discussed on the basis 
of the recent extension of the mode coupling theory for the liquid-glass transition to 
the model of the quenched-annealed binary mixture. It is in particular shown that, in 
confinement, the density correlation functions always decay to a non-zero infinite time 
value, even in the fluid state, and some clarification is given about the question of the 
relation between structure and dynamics in confined fluids. 
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In the past few years, there has been a strong interest for the dynamics of liquids 
under confinement and more specifically for their structural glass transition, which has 
been investigated by a variety of experimental techniques and by computer simulation 
[H Ej. One of the main goals of these studies was to improve our understanding of the 
concept of cooperativity, a key ingredient of many glass transition theories j^j. Indeed, 
there are now many evidences that, in the bulk, the dynamics of deeply supercooled 
liquids is strongly inhomogeneous and that correlated clusters of molecules, the so-called 
dynamical heterogeneities, play a crucial role in the slowing-down of the dynamics when 
the temperature is decreased. But, up to now, many aspects of the characterization of 
these dynamical heterogeneities have remained elusive. For instance, there is no clear 
consensus on their shape, their size and their evolutions with temperature. 

This is where the interest for confined glassforming systems came in. Indeed, 
confinement is a means to introduce geometrical constraints and new characteristic 
lengthscales (pore size, film thickness. . . ) in the system under study. Thus, by looking 
at the way the dynamics is modified under confinement compared to the bulk, one 
can hope to gain some insight on the properties of the dynamical heterogeneities. For 
instance, in the simplest scenarios, one expects from finite size effects a cutoff on the 
slowing-down of the dynamics as temperature is varied, when the typical size of the 
heterogeneities in the bulk would become larger than the characteristic lengthscale of the 
confining medium. Confinement would thus provide an indirect probe of the properties 
of the dynamical heterogeneities. 
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It turns out that the situation is more complex. Indeed, the previous line of 
reasoning requires that the physical phenomena which are specific of confined systems 
have a negligible impact on the dynamics of the imbibed fluid or at least that their 
influence is sufficiently well known that it can subtracted from the results. This is usually 
not the case and strong confinement effects are observed, in particular at the fluid-solid 
interface where structured layers of almost immobile molecules are often formed. So 
this is in fact the problem of dynamics under confinement as a whole which has to be 
addressed and not only at the level of a simple modification of the bulk dynamics. 

Theories dealing with the dynamics in confinement are far less developed than for 
the bulk and they usually are of a phenomenological nature [2]. There is yet a strong 
need for elaborate microscopic theories in this field. Indeed, the variety of systems 
to consider is immense. Porous media can differ in the shape, size, size distribution, 
connectivity, etc, of their pores. They can be made of various materials, leading to 
a wide range of fluid-solid interactions which adds to the already great variability 
of intermolecular interactions met with usual glassformers. A reasonable theoretical 
approach, able to catch many of these subtleties, would thus be of a great help. Applied 
to various models, it would allow to explore thoroughly the phenomenology of confined 
glassforming systems and maybe to disentangle the different physical effects which 
interplay in these systems. 

Recently, a step in this direction has been made with the extension of the mode 
coupling theory (MCT) for the liquid-glass transition jU El El El to a particular class of 
confined systems, the so-called "quenched-annealed" (QA) binary mixtures [Hj. In these 
systems, first introduced by Madden and Glandt P , the fluid molecules equilibrate in a 
matrix of particles frozen in a disordered configuration sampled from a given probability 
distribution. This class of models, which describe situations of statistically homogeneous 
and isotropic confinement, is thought to be able to reproduce most of the physics of fluids 
confined in materials like Vycor, controlled porous glasses or aerogels, and, in fact, the 
fluid dynamics and glass transition in some of its instances have been the subject of 
recent studies by molecular dynamics simulations [TU 1 ITf ) \T2\. 

In this paper, we will give a short presentation of the proposed extension of the 
MCT to QA systems (a more detailed account is given in reference jHj) and, on the basis 
of this approach, we will discuss a few aspects of the physics of confined fluids which 
are of relevance for the interpretation of experimental and computer simulation results. 

Before dealing with dynamics, one first has to consider some peculiarities of the 
statics of QA systems. QA mixtures are systems with quenched disorder, so that their 
theoretical description requires two types of averages, a thermal average denoted by 
(•••), taken for a given realization of the matrix, and a disorder average over the 
matrix realizations, denoted by TTT , to be taken after the thermal average. Like in 
the bulk, one is interested in the Fourier components of the microscopic fluid density, 
or, in short, density fluctuations, defined as p£(£) = Z)j=i e * qrj ^ > where q denotes 
the wavevector, Nf is the fluid particle number and Tj(t) is the position of fluid 
particle j at time t. A significant difference with the bulk is that, for a given matrix 
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realization, the translational invariance of the system is broken by the presence of 
the quenched component. This results in non-zero average density fluctuations at 
equilibrium, i.e., (p£) ^ 0. It is only after the disorder average that the symmetry 

is restored, leading to (pq) = 0, hence the description of the model as statistically 
homogeneous and isotropic. This property has a well known impact on the equations 
describing the structural correlations in such systems, for instance the so-called replica 
Ornstein-Zernike (OZ) equations [T3J Ej, where it leads to the splitting of the total 
and direct correlation functions of the fluid, h**(r) and c"(r), respectively, into two 
contributions, connected [h c {r) and c c (r)] and blocked or disconnected [h b {r) and c b (r)}. 
The separation of h**{r) into two terms leads to a similar property of the fluid structure 
factor SH = {pqP-^)/Nf = 1 + rifMf, where rif is the fluid number density, which can 
be expressed as S ff = S q + S b with S c q = ((p q — (p q ))(p^ q — (pLq)))/Nf = l + nfh c and 
Sq = (p<i)(P-q)/Nf = nfh b q , where f q denotes the Fourier transform of f(r). To fix all 
the notations, we define here the matrix-matrix and fluid-matrix structure factors and 
total correlation functions as well, which are given by 3™™ = (pqP™ q )/iV m = l + n m h q nm 

and S^ m = (p q p™ q ) / J NfN m = ^nfn m h^ m , where N m is the matrix particle number, 
n m is the matrix number density, and p™ = J2j=i e * qSj > where Sj is the fixed position 
of matrix particle j, is the q Fourier component of the quenched microscopic matrix 
density. 

The identification of two types of static fluid correlations has significant implications 
for the dynamics. Indeed, if one forgets for a moment the possibility of a dynamic al 
ergodicity breaking, one expects, using standard arguments, that lim^oo (pq(t)P-q) = 
(pq)(p_ q ), i.e., the normalized total density fluctuation autocorrelation function 4%(t) = 
(pq.{t)pL^) /(NfSH) does not decay to zero at long times, but rather 

lim <j> T (t) = -h. 

This general result only derives from the fact that the fluid is plunged in an 
inhomogeneous external potential, which here is due to the matrix. The existence 
of a non-zero limit thus is independent of the model and in particular it does not 
depend on the fluid-matrix interaction. To use the language of scattering experiments, 
it means that one should always expect to measure non- vanishing intermediate scattering 
functions or that an elastic contribution has to be present in the dynamical structure 
factors ITS] . We stress here that this is a true static phenomenon, and not one 
related to some ultra slow dynamical process. 

We can now turn to the dynamical theory. Following the preceding discussion, the 
proper dynamical variable to consider is the relaxing part of the fluid density fluctuations 
<5p q W = Pq(t) ~ (pQi rather than the full pl(t). Then, using standard projection 
operator methods [Sj, the equations of the MCT for QA systems are obtained [Hj. They 
are equations for the time evolution of the normalized connected density fluctuation 
autocorrelation function, 4> q (t) = (8 Pq(t)8 p_ q ) / '(NfSg) , and they consist of a generalized 
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Langevin equation, which is the same as for the bulk, 



+ {1% + & / drM g (t - r)0 9 (r) = 0, 



with Q 2 = q 2 kBT/(mSg), where m is the mass of the fluid particles, T the temperature 
and ks the Boltzmann constant, and of an expression for the memory kernel M q (t) = 

rv(t) + M g ( MC )(t), with 
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There are strong analogies between these equations and those for the bulk. Indeed, 
they share the same mathematical structure, so that most of their properties are already 
known from the extensive study of reference j3]. At a more detailed level, the expressions 



of the vertices V^. and V^, which determine the long time behaviour of <p q (t) and 
involve static quantities only, are very similar to those found in the bulk, for collective 
and tagged particle density relaxations, respectively 

About these vertices, a specific point is worth stressing. Indeed, one can see 
that the structural quantities on which they depend are not those related to the 
full density fluctuations, but the connected ones, which specifically characterize the 
static correlations of the relaxing part of the density fluctuations. This result, whose 
generality, we believe, goes beyond the model of QA mixtures, solves an apparent 
paradox mentioned in various places and sometimes considered as a challenge to a mode 
coupling description of the glass transition in confinement, that systems with identical 
structural properties can have significantly different dynamical behaviours [Tl"| I17j: In 
fact, only a nontrivial fraction of the total density fluctuations, a fraction which exists on 
top of the static average density fluctuations induced by the external potential in which 
the fluid particles evolve, is "active" in the determination of the dynamical properties of 
the confined fluid, and it is the evolution of this contribution which has to be considered 
to discuss the relation between the fluid structure and dynamics and not the total 
structure functions of the system. 

The solution of the MCT equations for QA mixtures offers no special difficulties 
compared to the bulk, since we essentially deal with the same type of equations. In 
reference [Ej, we have obtained by the methods of reference ^H] the dynamical phase 
diagram of a simple system consisting of a fluid of hard spheres confined in a matrix of 
hard spheres frozen in an equilibrium configuration, with both types of particles having 
diameter a j!2j . We present here typical examples of the time evolution of the predicted 
density correlation functions. Like in previous work, the system is characterized by 
two volume fractions <bt = 7rnta 3 /6 and <b m = nn. 



r(l) 



l cr 3 /6, and the necessary structural 




Figure 1. Connected density correlation function 4> q (t) for hard sphere fluids confined 
in matrices of identical hard spheres frozen in an equilibrium configuration, qa ~ 6.71 
is close to the position of the main maximum of the fluid structure factor. We use 
t mic = 160. (a) m = 0.2; from bottom to top: 0/ = 0.990}, 0.9990}, 0.99990}, 
0.999990}, 1.000010}, 1.00010}, 1.0010}, 1.010}, where 0} ~ 0.672 is the critical 
fluid compacity. The dotted line shows the asymptote of <fi q (t) at 0}. (b) m = 0.4; 
from bottom to top: 0/ = 0.90}, 0.990}, 0.9990}, 1.0010}, 1.010}, 1.10}, where 
0} ~ 0.279 is the critical compacity. 



quantities are computed using the Percus-Yevick approximation ^31 EI]- For simplicity, 
we neglect (j) q (t) in the generalized Langevin equation, so that the transient dynamics 
of 4> q {t) now has a time constant r q = t mic S q / (qa) 2 [T5] . 

Two types of ergodicity breaking transitions were found in reference jHj- For 
small matrix densities, a type B or discontinuous ideal glass transition line is met. It 
corresponds to the well known scenario found in the bulk, showing a two-step relaxation 
in the ergodic phase and a discontinuous change of the asymptotic value of (j) q (t) from 
zero to a finite value when the transition line is reached, corresponding to the divergence 
of the characteristic time of the second relaxation step. This is illustrated in figure 
la. For larger matrix densities, a type A or continuous ideal glass transition line is 
found. Here, one observes a one-step relaxation and the ergodicity breaking transition 
is continuous, i.e., the non-zero asymptotic value of 4> q (t) grows continuously from zero 
when the system enters in the non-ergodic domain. Such a behaviour is shown in figure 
lb. 

From the knowledge of 4> q (t), one can compute the total density correlation function 
4^{t), which is the quantity directly measurable in experiments and simulations. Both 
functions are indeed related by the simple formula 

m = §jut) + 4- 

Oq Oq 

In figure 2 are reported the results of this transformation on the data of figure 1. 
From these curves, difficulties can immediately be anticipated if one is to compare 
the predictions of the mode coupling theory with experimental or simulation data for 
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Figure 2. Total density correlation function (f^(t) for the same systems as in figure 
1. (a) <p rn = 0.2; <pfs as in figure la. The lowest dotted line shows the ratio S^/S^ 
at <fif, while the highest shows the corresponding asymptote of <f>q(t). (b) <p m = 0.4; 
4>fs as in figure lb. The dotted line shows the ratio S^/S^f at <j>^. 



<pq(t). Indeed, one sees that the blocked part of the static density correlations provides 
a density dependent background on top of which the glassy dynamics develops itself. It 
might thus be difficult to separate, in the long time behaviour of 4^(t), the evolutions 
which are of purely static origin from those which characterize the glassy dynamics of 
the system. This is especially critical in the case of type A transitions, since when the 
system enters in the ideal glassy state both contributions are, to first order, linear in 

<t>f-F f U 

In conclusion, we have demonstrated how the model of the QA mixture can be a 
very useful tool for a better understanding of the physics of confined glassforming liquids. 
Indeed, we have shown that this model, thanks to some crucial simplifying features, in 
particular the property of statistical homogeneity, allows for a rigorous illustration of 
some essential concepts which have to be taken into account when dealing with confined 
systems, like the distinction between static and relaxing density fluctuations. Now, with 
the corresponding extension of the MCT which is available, it seems that much progress 
could be done in our understanding of the dynamics of confined fluid systems. 
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